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O : Abstract 



I. INTRODUCTION 



At a future Linear Collider one will be able to determine the masses of charginos and neu- 
tralinos and their pair production cross sections to high accuracies. We show how systemati- 
*0 , cally including the cross sections in the analysis improves the measurement of the underlying 

mass parameters, including potential CP violating phases. In addition, we investigate how 
experimental statistical errors will affect the determination of these parameters. We present 
a first estimate on the lower limit of observable small phases and on the accuracy in deter- 
C\) ■ mining large phases. 

l> 

o 

On 

In the Minimal Supersymmetric Standard Model [[p the gauge boson and Higgs sectors are mapped to 
^ a set of neutral and charged fermions which mix to form neutralino and chargino mass eigenstates. The XX 
Q_i cross sections are parameterized by the gaugino masses M 1 , M 2 , the higgsino mass parameter fi, and the 
<—| ratio of vacuum expectation values tan (3. In CP non-conserving Minimal Supergravity models, the param- 
eter \i can be complex, and in the most general unconstrained Minimal Supersymmetric Standard Model 
the gaugino masses M 1; M 2 will also include phases; after rotating the wino fields, we can without a loss 
of generality choose 02 = 0, and are left with two additional phase parameters in the chargino/neutralino 
sector: 

/i^|/i|e^ Mi -> | Mi | e irpl (1) 

A future Linear Collider [Q] will enable the very successful LEP e + e~ precision analyses of the Standard 
Model to be extended so as to expose any underlying model. From the set of masses, total cross sections [Q] 
and distributions/asymmetries one can unambiguously determine the mass parameters for CP conserv- 
ing models. Yet it has not been quantified how well in a complete analysis one can measure small 
or large [[7p complex phases, and how this can improve their analytical determination from the physical 
masses Experimental bounds on electric dipole moments severely restrict the existence of phases 

which do not rely on large cancellations, allowing only small variations due to non-gaugino MSSM param- 
eters, like the trilinear couplings A+ [|],^p. Phases of the order O(ir/I0) or even smaller, like 0(7r/lOO), 
hardly influence naive CP conserving observables like masses and total cross sections and will be diffi- 
cult to determine at a Linear Collider. We will give a first estimate of when statistical errors render small 
phases unobservable and how well large phases can be measured, taking into account realistic experimental 
uncertainties. 
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II. ANALYSIS OF MASSES AND CROSS SECTIONS 



To estimate the effect of CP violating phases on masses and total cross sections, we examine five 
scenarios. Three of them are derived from a set of mSUGRA parameters taking m = 100 GeV and 
m 1 /2 = 200 GeV [0, adding small phases (pi and Since the phase (fi 1 = O.Oln in scenario (2) is very 
small, this parameter set is similar to an mSUGRA scenario. The large phase scenario [0] avoids constraints 
from electromagnetic dipole moments through cancellations. [] In the wino lightest supersymmetric 
particle model pTTp, additional phases are attached to the largest entries Mi,// in the neutralino/chargino 
mass matrices. Simply adding small phases to CP conserving models will reveal how well we will be able 
to distinguish CP conserving and CP violating models. Only in the large phase scenario (4) are the phases 
a basic feature of the model: 
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Mismeasuring masses and total cross sections within their experimental uncertainties is the limiting 
feature expected to spoil the determination or even observation of phases, because the numerical effect, 
especially of small phases, is generically small. For each of the five scenarios we calculate the set of 
chargino/neutralino masses and the pair production cross sections for all possible final states. We regard 
this sample as a set of experimental observables, mismeasured with a given Gaussian uncertainty. Using 
either a global fit or the inversion algorithm described in the Appendix, we then extract again the parameters 
Mi, M 2 , //, tan (3 and <f>i, M from the smeared set of observables. 

For the design parameters of the future Linear Collider we choose two possible levels of performance: 
a 0.5 TeV machine collecting an integrated luminosity of 0.5 ab _1 (500 fb _1 ), and a high performance 
version with 1 TeV and 1 ab _1 . Although the smaller machine might not be able to produce large numbers 
of higgsino pairs, the high performance version will extend its mass reach only at the expense of the gaugino 
cross sections which drop like 1/s, as illustrated in Figure [I]. In Table | we present the set of masses and 
cross sections for the scenarios (1) and (4). 

Typical errors for the mass measurement at a 800 GeV Linear Collider have been estimated for a scenario 
very similar to scenario (1), rendering 0.05/0.07/0.3/0.6 GeV for the four neutralinos and 0.035/0.25 GeV 
for the charginos [Q|. We use these one sigma error bars for all scenarios, because the typical mass scales 
are similar. Moreover, we disregard additional theoretical errors or higher order corrections; the latter 
merely shift the theoretical predictions and affect the currently unknown errors. The effect of mismeasured 



t-channel slepton masses on non-polarized total cross sections is negligible [|12|]; we assume their exact 
determination from the direct production to reduce the number of observables varied in the analysis. The 
total cross section for x\ P arr production is not part of our set of observables. 



1 The very small value of tan (3 is neither a generic feature of the model [10] nor of our analysis. However, we prefer 
to quote the exact values of Ref. [[7p. 



2 



For the fits, we minimize x 2 = Si(^reconstr,i — ^mcas.i) 2 /^ using MINUIT. We choose a Gaussian 
probability distribution for the observed values for masses and cross sections x meaSj i, disregarding systematic 
errors. Assuming an efficiency e = 10% we obtain e\ = a-J(eC) for the cross section measurements. The 
result of the inversion is a set of reconstructed theoretical input parameters x reconstr; i. Randomly varying all 
observables gives 10000 best fitting sets of reconstructed input parameters, displayed in each plot. As long 
as the result is peaked we give the standard deviation RMS of the pseudo-measurements of the two phases. 

Modified mSUGRA: Gaugino mass unification M 2 ~ 2M± together with radiative electroweak symme- 
try breaking fi ^> Mi, M 2 yields a mass hierarchy in which the heavy higgsinos have generically smaller 
cross sections than the light gauginos; their largest cross section is < 200 fb for xf pair production. Naively 
one might guess that in particular the neutral higgsinos would have large cross sections, since they couple 
to the s-channel Z boson, whereas the gauginos are produced through t-channel slepton exchange with the 
sleptons being heavier than the Z. In fact, the mixing and the phase space factors dominate the hierarchy 
of the cross sections. Thus we have to measure small cross sections of ^ 70 fb to determine M . Moreover, 
large errors on the higgsino masses do not allow as precise a determination of 0^ as of <f>\. Figure ^ shows 
that 0i = 0.17T can be determined to ~ 30% with the 0.5 TeV machine, whereas the best fit values for 
M = 0.17T render a RMS value of 46%. For the high performance design the M determination becomes es- 
pecially difficult. This cannot be improved by applying cuts on % 2 ; all 10000 best fits are of similar quality. 
The x 2 value in this naive approach contains only little information on how closely the reconstructed values 
lie to the theoretical input ones. 

A very small phase X < 0.017T will be indistinguishable from zero: the distribution of the fits in Figure ^| 
does not allow one to extract a non-zero central value of <f>\. Moreover, the 0^ distribution develops a second 
peak for zero phases. This indicates that there are two possible minima in x 2 , leading to best fits of similar 
quality. Since this is a drawback of the weak dependence of masses and cross sections on small phases, the 
fits of the mass values in Figure [| still vary by only (9(lGeV)Q A small value for tan (3 can be determined 
to ~ 10% for both of the collider performances. 

The case of tan (3 = 30 in Figure |] shows a similar behavior to Figure 0. There is, however, a large 
uncertainty in tan/3; the steep rising of the tangent is not reflected in the sine and cosine behavior of the 
mass matrices and yields a relative error of > 20% on the measurement of large tan (3 values. This tan (3 
uncertainty influences the measurement of both of the phases, but X = O.ln = 0^ can still be distinguished 
from zero and probably be measured at a 0.5 TeV collider. 

For all mSUGRA type scenarios, the larger number of observed processes at the high energy collider 
does not increase the accuracy of the measurements. As depicted in Figure [I], gaugino pairs have large cross 
sections at the 0.5 TeV machine, and the phase of the higgsino parameter has a sufficient impact on the 
X1X2 production cross section and the higgsino masses. Even compensating for the smaller cross sections 
by doubling the luminosity of the 1 TeV collider does not improve the results of our fits. The additional 
observables are not sensitive enough to the phases to effectively contribute to the quality of the fit. However, 
this only holds as long as the threshold for X1X2 production is slightly below 0.5 TeV. 

Large phases: Avoiding the dipole moment constraints by fine tuning ji leads to very heavy higgsinos in 
the given scenario (4). They will mainly be produced together with a gaugino. The set of observable cross 
sections ranges up to ~ 170 fb, even for the high energy design. Large phases 1; 0^ lead to a considerable 



2 We have also fitted a zero phase mSUGRA scenario 0i = = M leaving the phase values free. The obtained 
RMS values are somewhat dependent on the allowed frame in the fit; typical RMS values are 0.03/0.04 for the 0i 
distribution and 0.06/0.07 for M assuming the low/high performance collider design. 
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effect on masses as well as on cross sections. The phase 4>\ can be measured to ~ 5(7)% at the low (high) 
performance machine, see Figure |5|. The lower energy removes the higgsino masses and cross sections from 
the sample; this renders the parameter fi indeterminable. The distribution peaks around a wrong central 
value, again showing that there is no sensitivity to the phase. On the other hand the accuracy of the Mi and 
0i measurement is not affected. Using 1 ab _1 at a 1 TeV collider improves the measurement to % 5%. 

Since fitting to the large phase scenario yields the best measurements of the phases we also try to extract 
the theoretical parameters from the masses alone. In the right column of Figure |5] we compare the fit to 
the masses to the result of the algorithm described in the Appendix. Neither the fit to the masses nor the 
algorithm lead to a measurement of (f>i\ the distribution of values extracted from the algorithm is entirely 
flat. The corresponding fit would prefer arbitrary wrong minima, rendering the result completely dependent 
on the setup of the fitting procedure. Both approaches yield a similar central value and width for only 
because the fake minima for the complete set of parameters give identical M values. This example shows 
that only a large number of observables can guarantee a reliable determination of phases, compensating for 
the generally weak dependence of the observables on the phase parameters. 

Modified wino LSP: Potentially dangerous features of the wino LSP scenario are large values for Mi 
and fi, the complex entries in the mass matrix. The phases could effectively decouple. This, however, 
merely affects the accuracy of the fitted Mi values displayed in Figure ^. It is limited by the large errors of 
the heavy neutralino/chargino mass measurements. A small absolute value of /x on the contrary can easily 
be determined from the masses alone. Due to strong mixing in the mass matrices, <pi can be measured with 
% 20 percent uncertainty, whereas fitting the phase relies on the small mixed cross sections and might 
only distinguish a finite value from zero. 

IH. SUMMARY 

For several different scenarios we estimate a possible straightforward measurement of CP violating 
supersymmetric phases. Small phases are added to mSUGRA models and to a wino LSP model: typical 
values around n/10 for 1; M can be determined from total cross sections and masses by simple fits; they 
will not be hidden by anticipated statistical errors. However, very small phases, of the order of 7r/100, as 
preferred by electric dipole moment analyses [0,^], can hardly be distinguished from zero. In this case, the 
sample of observables behaves like a CP conserving set, and the mass parameters and tan [3 are accessible to 
a high degree of accuracy [f|]. Large phases can easily be determined if the corresponding gaugino/higgsino 
states are produced. For the model under consideration [|7]] the set of masses and cross sections at a 1 TeV 
collider can be used to determine the values to a few percent. However, the set of masses alone is heavily 
affected by experimental errors. There still is a remaining sign ambiguity of the phases. Forward-backward 
asymmetries, left-right asymmetries, and other explicitly CP non-conserving observables will considerably 
improve our estimates A more sophisticated analysis including some detector simulation would have 
to be carried out to determine where the lower limits for detectable very small phases will finally lie. 
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APPENDIX A: INVERSION ALGORITHM 



From the complete set of neutralino/chargino masses we can unambiguously determine the underlying 
mass parameters.^ In the following we denote the absolute values \M\\, |/i| by Mi, fi. First we calculate all 
parameters as a function of M 2 : 
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Here F = My / /{2M 2 fj), Fi = M 2 (? w FqI (Mis 2 ), and s w , c w stand for the sine/cosine of the weak mixing 
angle. After computing these parameters, 

u C l + F l sm 2 2(3 

COS 0i + M = — 

sm 2(3 

serves as a self consistency relation and thereby determines M 2 . 



3 Ways have been shown to extract these parameter from a smaller set of masses |9||. However, the complete set is 
well suited to illustrate the result from the fits. 
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TABLE I. Neutralino/chargino masses and cross sections at a 1 TeV Linear Collider for the modified mSUGRA 
scenario (1) and the large phase model (4). The large X1X1 cross section is not used in the analysis. 
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Figure 1. The largest cross sections for the mSUGRA type parameters (1) and <p\ = 0.17T, M = 0.17T as a 
function of the collider energy. The solid lines indicate cases where CP violating phases do not influence the explicitly 
CP conserving observable. 
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Figure 2. Parameters obtained from best fits to masses and cross sections with the mSUGRA type parameter 
choice (1) and <pi = O.Itt, (p^ = O.Itt. All central values are indicated by arrows on the axes. The grey (green) and 
black points correspond to the low and high performance collider designs, respectively. The straight boundary in the 
M2-tan (3 plane is an artifact from restricting the phases to [0, it]; this makes it easier for MINUIT to find a global 
minimum in x 2 - 
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Figure 3. Parameters obtained from best fits to masses and cross sections with the mSUGRA type parameter 
choice (2) and = 0.0l7r, m = 0.17T. The central values are indicated by arrows on the axes. The grey (green) and 
black points correspond to the low and high performance collider designs, respectively. 
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Figure 4. Parameters obtained from best fits to masses and cross sections at the large tan j3 mSUGRA parameter 
point (3) with <pi = O.ln, = O.Itt. All central values are indicated by arrows on the axes. The grey (green) and 
black points correspond to the low and high performance collider designs, respectively. 



8 



(0.5TeV, 0.5atrV(lTeV, lab" 1 ) 



6 masses alone: algorithm/fit 



78 

76 

74 

72 
480 

460 

440 

420 



MJGeV] 


1 










~ , , , , , , t , 





0.4 



0.5 



0.6 



|J.[GeV] 

'.VV V .• 




0.7 



-RMS=0.023/0.024 



0.8 



0.9 



2000 



10000 best fits 



>,/j[ 



0.4 



0.5 



0.6 



2000 




80 



78 



76 



74 



72 



451 



450 



449 



448 



500 



M.[GeV] 




>,/j[ 



0.5 



0.75 



_ HfGeV] 



2000 



0.5 



mean=0.72/0.79 
RMS=0.09O/0.O75 



0.75 








' , , , t , , , 


<t>/t 






0.7 0.8 0.9 0.5 0.75 1 

Figure 5. Parameters obtained from best fits to masses and cross sections (left) and to masses (right) at the large 
phase parameter point (4) 4>\ = 0.57T, 0^ = 0.87T. The central values are again indicated by arrows on the axes. The 
grey (green) and black points correspond to the low and high performance collider designs, respectively. For the grey 
(blue) dots in the right column of plots the fit to the masses is replaced by the algorithm described in the Appendix. 
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Figure 6. Parameters obtained from best fits to masses and cross sections at the wino LSP parameter point (5) 
with 4>\ = 0.17T, (ftp = O.Itt. The central values are indicated by arrows on the axes. The grey (green) and black points 
correspond to the low and high performance collider designs, respectively. 



10 



